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ABSTRACT 

We present rest-frame NIR luminosities and stellar masses for a large and uniformly-selected pop¬ 
ulation of GRB host galaxies using deep Spitzer Space Telescope imaging of 119 targets from the 
Swift GRB Host Galaxy Legacy Survey spanning 0.03 < z < 6.3, and determine the effects of galaxy 
evolution and chemical enrichment on the mass distribution of the GRB host population across cos¬ 
mic history. We find a rapid increase in the characteristic NIR host luminosity between 2 ~ 0.5 
and z ~ 1.5, but little variation between z ~ 1.5 and z ~ 5. Dust-obscured GRBs dominate the 
massive host population but are only rarely seen associated with low-mass hosts, indicating that mas¬ 
sive star-forming galaxies are universally and (to some extent) homogeneously dusty at high-redshift 
while low-mass star-forming galaxies retain little dust in their ISM. Comparing our luminosity dis¬ 
tributions to field surveys and measurements of the high-z mass-metallicity relation, our results have 
good consistency with a model in which the GRB rate per unit star-format ion is constant in galax¬ 
ies with gas-phase metallicity below approximately the Solar value but heavily suppressed in more 
metal-rich environments. This model also naturally explains the previously-reported “excess” in the 
GRB rate beyond z > 2; metals stifle GRB production in most galaxies at z < 1.5 but have only 
minor impact at higher redshifts. The metallicity threshold we infer is much higher than predicted 
by single-star models and favors a binary progenitor. Our observations also constrain the fraction of 
cosmic star-formation in low-mass galaxies undetectable to Spitzer to be small at z < 4. 

Subject headings: gamma-ray burst: general — galaxies: star-formation — galaxies: evolution 
galaxies: photometry — galaxies: high-redshift 


1. INTRODUCTION 

Long-duration gamma-ray bursts provide a unique 
and powerful means of studying galaxy evolution: as 
the extremely luminous explosions of massive stars 
ijWooslev fe Blooii] [20061) . they are d etectab le out to 
very high redshifts dTanvir et al.ll2009t ISalvaterra et al.1 
l2009l ICucchiara et all 1201 lh . provide backlights with 
which to st udy the intervening ISM a nd IGM in 
detail (e.g., iProchaska et al.112007L 120091 : iFvnbo et al.1 
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l2009f) and—perhaps most importantly—their rate is 
governed (at least in part) by that of star-formation. 
The redshift distribution of GRBs and the proper¬ 
ties of their hosts can therefore be used to study 
the evolution of the cosmic star- f ormat ion rate den¬ 
sity wit h ti me (e.g., »l ot anil ll997l IWiiers et al.l 119981 : 
iBlain fe Nataraiaiil l200Ct iPorciani fe Madaul 120011) and 

characterize the sites where stars were formed in the 
early universe, independent of many of the foils that hin¬ 
der traditional techniques such as uncertain dust cor- 
rcctio ns or incompleteness t o low-luminosity galaxies 


(e.g., Diorgovski et al.l 2001c Ramirez-Ruiz et al.1 

2002 

Tanvir et al.l 20121: 

Trenti et al.l2012i: Schulze et aU 

2015 

Greiner et al.l 2015 

)• 


However, only a tiny minority of massive stars actu¬ 
ally produce a gamma-ray burst, and as a result the 
GRB rate may in principle be influenced significantly 
by factors other than purely the rate of star-formation, 
such as metallicity. A detailed understanding of any 
such variations is essential to apply GRB-based infer¬ 
ences of the cosmic star-formation rate density or other 
broader topics with confidence. Furthermore, they would 
serve to illuminate the nature of the GRB progeni¬ 
tor: the simplest single-star progenitor formation sce¬ 
narios imply a strong preference or even requirement for 
very low metallicities, due to the need to avoid loss of 
mass and angular momentum in line-d ri ven winds (e.g., 
Heeer et al.l l200d iHirschi et al.l 120051: lYoon fe Langerl 

2005HYoon et al.ll2006l : iWooslev fe Hegerll2006l) . Binarv- 
evolution models are more flexible and typically im- 
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ply relatively modest metal sensitivity le.g.Jlzzard et all 
1 2004 iFrver fc Hegedl2005l : iPodsiadlowski et al.ll2010l h 


The extent to which the GRB rate varies as a func¬ 
tion of environmental factors such as metallicity remains 
an open question observationally. Studies of the red- 
shift distribution alone provide some insight: the co¬ 
moving GRB rate density exhibits qualitatively simi¬ 
lar behavior as the star-formation history over most of 
cosmic time, but shows a modest surplus of GRBs at 
high-redshift (z > 2 — 3) relative to what would be 
expected if the GRB rate tracked the star-formation 
rate exactly (e.g., iKistler et al l 1200 81: iBut.le r et alJ 2£)_Hjj; 
Wanderman fc Piranl l201fll: iRobertson fc Ellisl 20T2I: 


Jakobsson et al. 120121 iLien et al.l I20l4 iPerlev et alJ 
2016f) . This alone seems to rule out extremely metal- 


dependent production models and provides some evi¬ 
dence in favor of a more modest dependence. However, 
this technique is necessarily imprecise: the comoving 
GRB rate measured at any redshift is an average over 
all galaxies at that epoch (a diverse population span¬ 
ning orders of magnitude in metallicities, star-formation 
rates, etc.) and cannot reveal which types of galaxies 
contribute most (or least) to the GRB rate. A thorough 
understanding of the link between GRB production and 
environment requires also characterizing of the popula¬ 
tion of GRB host galaxies directly. 

Many host-galaxy studies have been performed over 
th e past decade; a su mmary of these efforts can be found 
in IPerlev et al.l 120161 (h ereafter. Paper I) and in the re¬ 


view of lLevesqcui 1201 l i. Nearby (z < 1) GRBs are com¬ 


monly found in hosts that are blue, young, low-mass and 
metal-poor in comparison to the star-forming field pop¬ 
ulation, but not necessa ril y extr eme in regards to any of 
these attributes dSavaglio et al.l 1200911 . They occur only 
very rarely in metal-rich and massive galaxies in this red- 
shift range: large, massive spirals with super-Solar metal¬ 
licity similar to the Milky Way contribute substantially 
to the local star-formation rate d ensity but are almost 
never observed to host GRBs (e.g.JLe Floc’h et al.ll2006t 
iGraham &; Fruchterll20 1 3l : iVergani et al.ll2015lh 

The situation at higher redshifts is much less clear. 
If the rate of GRBs is controlled by metallicity, the de¬ 
clining metallicities of galaxies with increasing redshift 
should result in a high-redshift host population that 
is much more repre sentative of star-formation overall 
comp ared to z ~ 0 dFvnbo et al.ll2006t iKocevski et alJ 
I2009T) . The first deep NIR and Spitzer studies of GRB 
hosts at z > 1 produced few detections at these wave¬ 
lengths, suggesting a host population that remained 
quite low in average mas s (Lc Floc’h ct al. 2003, 120061 : 
iCastro Ceron et ahf 120101 : lLaskar et alJ 1201 111 . How¬ 
ever, these early samples generally under-represented 
or omitted dust-obscured bursts: more recently, signif¬ 
icant numbers of GRBs in massive and metal-rich galax¬ 
ies have been fou nd by targeting dust-obscured “dark” 


bursts specifically (iKriihler et al. [gOPi iRossi et al.ll2012t 
IPerlev et al.ll2018fc IHunt et al.ll2014li 41so. even without 
specifically including dark bursts, the UV luminosity dis¬ 
tribution of small samples of GRB hosts at z = 2 — 4 
seems to show reasonable consistency with a popula- 


(iFvnbo et al. 

20Q^; Chen et al. 20091; Greiner et al. 2015; 

Schulze et al 

20H). 


Actually extending these results to firm conclusions 


about the properties of the high-z host population or 
the nature of the progenitor is challenging: high-redshift 
massive GRB hosts clearly exist, but their actual abun¬ 
dance, and where they fit in the overall distribution 
of properties of host galaxies at z > 1, is diffi cult to 
quantify, and the t opic is still debat ed (Savagliol 120121 : 
IPerlev et alJ 120131 : IHunt et all 120141 . Observational 
studies conducted so far have been heterogeneously- 
targeted (biased for or against obscured bursts, and 
sometimes biased in favor of bright and luminous hosts 
that are easiest to identify and study), limited in depth or 
wavelength coverage (based on observations in only one 
or two filters and difficult to connect with real, physi¬ 
cal host properties), and/or too small or too limited in 
redshift to characterize the high-redshift population. 

In Paper I we introduced a new survey of the GRB 
host galaxy population (the Swift Host Galaxy Legacy 
Survey, or “SHOALS”), a new project designed to move 
beyond these limitations. Our survey provides the first 
host galaxy sample that is unbiased (homogeneously tar¬ 
geted) and sufficiently large (119 objects) to statistically 
examine redshift evolution in the host population in de¬ 
tail, and we are actively observing the sample at many 
different wavelengths and to sufficient depths to detect 
and thoroughly characterize the properties (M*, SFR, 
etc.) of each individual host. 

In this paper we present our Spitzer 3.6 pm. photom¬ 
etry of the SHOALS sample. Spitzer observations are 
particularly key to the effort, since the rest-frame NIR 
luminosity probed by Spitzer directly pr obes the stellar 
mass—which in turn is tied (e.g., iTremonti et al J l2001ll 
to the host metallicity. Metallicity has, for both theoret¬ 
ical and observational reasons, traditionally been viewed 
as likely to be the primary driver controlling the GRB 
rate as a fraction of the rate of star formation (the GRB 
efficiency), although it is not the only candidate. 

Our Spitzer-IRAC observations and data analysis are 
described in fJ3J supported by the redshift measurements 
and higher-resolution ground and space-based imaging 
described in Paper I. We outline our results in fJ31 show¬ 
ing the near-IR luminosity distribution (a good proxy 
for the mass distribution) of GRB hosts and its evolu¬ 
tion with redshift from the local Universe out to redshift 
6. We discuss the connection between host properties 
and the degree of obscuration of the afterglow, and its 
implications for future surveys as well as the nature of 
the ISM in high-redshift galaxies. We attempt to model 
the redshift evolution of the host population using a sim¬ 
ple luminosity-dependent cutoff model and show that the 
near-IR properties and redshift distribution of the pop¬ 
ulation are well-described by a simple model in which 
GRBs are heavily suppressed above a maximum metal¬ 
licity (oxygen abundance) t hreshold of 12_+ log[0/Hl = 
8.94 ± 0.04 using the scale of lKobulnickv fc Kewlevll2004 
(KK04). We compare our results to previous work and 
discuss implications in |J4] 

2. OBSERVATIONS AND ANALYSIS 
2.1. IRAC Observations 

The Spitzer Infrared Array Camera (I RAC. Ifazio et all 
l2004fl is a four-channel, mid-i nfrared imager opera ting on 
the Spitzer Space Telescope (jWerner et al.l 1200 4). Since 
the end of the cryogenic mission in 2009, only two of 
the channels are functional: channel 1, operating be- 
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Fig. 1. — Spitzer-IRAC 3.6/xm imaging of all 119 uniformly-selected GRB host galaxies from the SHOALS sample. The red circle denotes 
the best-available position of the GRB afterglow from the sources identified in Table 3 in Paper I; the yellow position is centered on the 
GRB host galaxy (Table [2] the center location is measured from optical imaging where possible, otherwise from the Spitzer data itself) 
or on the afterglow location if no host galaxy is detected at any waveband (the latter case is represented by a dashed aperture). Due to 
Spitzer’s depth and large PSF, many of the host galaxies are moderately to severely blended with neighboring objects. Images are 9 " 
across; North is up and East is left. 
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2.— Same as Figure[T] but after subtraction of nearby sources using our galf it-based scripting procedure to isolate the host galaxy. 
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tween 3.2-3.9 ^m; and channel 2, operating between 4.0- 
5.0/rm. (In this paper we will refer interchangeably to 
the channels using their channel numbers and central 
wavelengths: 3.6 fj,m for channel 1 and 4.5 /im for channel 
2 -) 

Prior to the start of our project, we examined the 
Spitzer observations catalog to determine which events 
within the SHOALS uniform sample (outlined in detail 
in Paper I) had previously been observed by the tele¬ 
scope. A large number (almost half) had observations 
already present in the archive (see Table [2] for details). 
We did not request re-observations of these targets. For 
the remaining targets we conducted new observations as 
part of our Cycle 9 Large Program (GO-90062; PI Per- 
ley)- 

Our observational strategy was based on the redshift as 
measured at the time of the start of our campaign (Fall 
2012), with deeper observations used for more distant 
targets. Integration times were typically 0.2-0.5 hours 
for z < 1; 1.5 hours for 1 < z < 2; 3 hours for 2 < 
z < 3, and 6 hours for z > 3. For events at unknown 
redshift at the time of the proposal, we integrated for 
1.5 hours. In each case a 100 s frame time was used, 
dithering using a medium-step dither pattern except in 
a few cases where the desire to avoid a nearby bright 
star favored a more compact dithering arrangement. For 
targets at unknown redshifts and for those thought to 
be dark/dust-obscured, we also obtained observations in 
channel 2 (4.5 fi m) using the same exposure prescription 
as for the channel 1 observations. 

We downloaded the PBCD (Post-Basic Calibrated 
Data) observations as they became available in the 
Spitzer Legacy Archive. We use the default astrometry 
provided by the pipeline without further additional align¬ 
ment, effectively establishing the IRAC images (which 
are aligned against 2MASS with an accuracy of 0.3" by 
the IPAC pipeline) as the reference system for the survey. 


2.2. Deblending 

Source confusion and contamination are common in 
even moderately deep Spitzer imaging due to the instru¬ 
ment’s large PSF (~ 1.8"). Indeed, inspection of our im¬ 
ages (Figure[l]) demonstrates that a significant fraction of 
the target host galaxies are contaminated to some degree 
by flux from neighboring, unrelated objects. Accurate 
photometry requires removing this contamination. 

We manually examined each Spitzer image simultane¬ 
ously with any higher-resolution ground-based or HST 
imaging available of the field (see Paper I) to identify 
the host-galaxy target and any nearby objects which may 
contaminate the host or sky apertures. Using a custom 
DS9/python script, we drew a fitting box around each 
host position and estimated the locations and sizes of all 
detected sources inside it (as well as those of any bright 
sources outside the box edges that contribute significant 
flux within the box). These were then sent as i nitial in¬ 
puts via a python wrapper to galfit dPeng et all 12002 ! 
to iteratively fit and subtract all sources within the box. 
The default PBCD image is used as the fitter input im¬ 
age with the PBCD uncertainty map providing the pixel 
uncertainties. For the input PSF we use our own merged 
combination of the core and extended PRF files from the 


IRAC website 13 , with the final PSF oversampled by a 
factor of 5 relative to the PBCD image resolution (or 
by a factor of 10 relative to the native pixel scale.) We 
use a Sersic model with index fixed at n = 1 (typically 
circular, but elliptical if the source is visibly elongated) 
or a point-source model depending on the appearance 
and brightness of the target. The residual images were 
inspected after the fit, and the routine was re-run with 
new inputs if necessary until all sources near the host 
galaxy (including the host itself, if well-detected) were 
removed as completely as possible from the final residual 
map. We then created a final output image by construct¬ 
ing a galfit model using the best-fit parameters for all 
components except for those associated with fitting the 
host itself, and subtracted this model from the initial 
PBCD image to isolate the host galaxy. The subtracted 
maps are shown in Figure [5] This procedure was effec¬ 
tive at removing contaminating flux for all of our fields, 
usually leaving negligible residuals except in the cores of 
very bright stars and galaxies. 

2.3. Host-Galaxy Photometry 

Our galfit procedure returns calibrated PSF-based 
magnitudes using the appropriate zeropoints for each im¬ 
age. For our final photometry, however, we employ the 
procedure recommended in the IRAC handbook of calcu¬ 
lating aperture photometric magnitudes. These are cal¬ 
culated via the subtracted images described above using 
a custom IDL wrapper around the aper procedure in the 
Astronomy User’s Library, employing an aperture radius 
of 1.8" (1.5 native pixels or 3 resampled pixels). Zero- 
points (including aperture corrections) are based on the 
corresponding values reported in the Spitzer-IRAC hand¬ 
book, interpolated using our PSF model 1 N2.2I) to permit 
the use of fractional native-pixel aperture sizes not given 
in the handbook. Apertures are centered on the host- 
galaxy localization as determined from our ground-based 
imaging or from the fitting/subtraction procedure; the 
aperture centers employed are given in Table [2] We esti¬ 
mate uncertainties via the scatter of pixel values in the 
sky annulus (within the aper procedure), rescaled using 
an empirical correction for correlated noise expected in 
subpixel-sampled images. We measured this correction 
to be a factor of 1.25 (in flux) for our PBCD images by 
placing a large number of random apertures in blank re¬ 
gions of several non-confusion-limited images and calcu¬ 
lating the true scatter in the resulting flux values. Error 
estimates do not explicitly include systematics associated 
with the subtraction procedure, which are expected to be 
minor. 

In many cases the host galaxy is not detected in the 
IRAC imaging. For these cases, we calculate an upper 
limit. If the host is detected in ground-based imaging we 
fix the aperture location at the host position and calcu¬ 
late a 2 <t upper limit on the flux at that location. If we 
do not know the host location, we center the aperture on 
the best-available afterglow position and calculate a 3 a 
upper limit. (The higher er threshold reflects the addi¬ 
tional uncertainty associated with not knowing exactly 
where to place the aperture.) 

Photometry on some of the hosts presented here 
were also published in previous work (in particular 

13 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/calibrationfiles/psff 
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Fig. 3. — Near-infrared absolute magnitudes of a uniform sample of 119 GRB host galaxies (110, or 93%, with measured redshift) 
and field galaxies. Magnitudes are AB and correspond to a rest-frame wavelength of A re st = 3.6 ^im/(l-\-z), the rest-frame equivalent of 
the IRAC channel 1 filter. Green points indicate bursts not known to be obscured; red points indicate dust-obscured and “dark” bursts. 
Darker shades of both colors indicate redshifts measured from late-time host galaxy observatio ns, while redshif t s of l ighter-shaded points 
are measured from target-of-opportunity afterglow observations. /C s -selected field galaxies from fKaiisawa et alJ (|2011f) are plotted in gray 
with area scaled according to SFR. The right panel shows GRBs at unknown redshift (arbitrary x-axis); in this case colors indicate the 
source of the redshift upper limit. Most of these limits are close to m ^ 25, which is translated to the main panel as an orange curve. 
GRB 070808 is shown at z = 1.35. ’’The cyan curve shows the redshift-dependent median magnitude of the GRB host sample, with the 
shaded region denoting the la bootstrap uncertainty on this value. 


lLaskar et alJ 1201 ll and iPerlev et alJ I2013TI but in all 
cases the analysis here is independent, using the semi- 
automated subtraction and photometry procedure above 
consistently for all targets. Our photometry is gener¬ 
ally consistent with these previously published values. 
We note that further observations and analysis (better 
de-confusion and precise localization of additional host 
locations using deeper images) may provide additional 
improvements in the future so the values reported here 
do not necessarily represent the final values for the sur¬ 
vey, but we do not expect large changes and the current 
photometry should be adequate for the purposes of this 
paper. 

2.4. Mass Conversion 

The luminosity of a galaxy as measured in the IRAC 
3.6/im band is a good tracer of stellar mass even out 
to high redshifts, since it probes emission redward of 
the Balmer break (dominated by evolved stars and only 
moderately dependent on population age) at nearly all 
redshifts. While a precise stellar mass estimate requires 
an SED fit to many photometric data points (to incorpo¬ 
rate age, extinction, and similar dependencies), a reason¬ 
able estimate can be obtained from Spitzer observations 


alone. 

For each of our (known-redshift) galaxies, we calcu¬ 
late the the quantity M 3 . 6 /(i +z ) = m 0 bs, 3.6 - n{z) + 
2.51og 10 (l + z); this is the luminosity (AB absolute mag¬ 
nitude) of the host galaxy at the rest wavelength ob¬ 
served by the 3.6 gm filter after cosmological redshift, 
i.e., the absolute magnitude at a wavelength of 3.6/(1+ z) 
fim. This quantity alone can be used as a reasonable 
stellar-mass proxy, but we also calculate stellar mass esti¬ 
mates directly using the following prescription. For each 
decade in stellar mass (1O 7 M 0 , 1O 8 M 0 , etc.) and for a 
grid of different redshifts between z =_0 — 10 . we cre- 
ate a model galaxy SED by summing iBruzual fc Charlotl 
(120031) galaxy SED templates for a continuous star- 
formation history between the present time and a maxi¬ 
mum age determined by the mass (100 Myr for < 1O 8 M 0 , 
300 Myr for 1O 9 M 0 , 700 Myr for 1O 1O M 0 , and 3 Gyr for 
> lO 11 Af 0 ; if this age is greater than the age of the Uni¬ 
verse at the template redshift we instead use the age of 
the Universe as the maximum age). We assume mod¬ 
est dust attenuation {Ay of 0.05 mag for < 1O 8 M 0 , 0.1 
mag for 1O 9_1 °M 0 , and 0.3 mag for > 10 n M 0 ) 14 and a 

14 These preferred values of Ay and age were chosen by compar- 
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Fig. 4.— Cumulative luminosity distribution of GRB hosts in a 
variety of redshift bins spanning most of cosmic history. Dashed 
lines indicate upper limits (upper limits are plotted as 1 magnitude 
fainter than the measured limiting value; this portion of the curve 
should not be taken as representative of the true distribution.) The 
host galaxy luminosity function increases sharply from z ~ 0.25 to 
z 1.5 and then changes relatively little out to higher redshifts. 

IChabrieil (l2003f ) IMF. For each SED we then calculate via 
synthetic photometry the same quantity which we mea¬ 
sure using Spitzer, M 3 .6/(i+ z ). The observed magnitudes 
can then be translated to stellar masses via interpolation 
on the mass-redshift grid. 

Compared to, e.g., scaling based on the rest-frame K- 
band luminosity, this empirical procedure avoids the use 
of fc-corrections and allows for the expected not-fully- 
linear scaling of luminosity with mass (because lower- 
mass galaxies tend to be younger and have lower mass- 
to-light ratios.) Of course, like all mass derivation proce¬ 
dures, it is subject to various assumptions related to the 
choice of stellar templates, IMF, star-formation history, 
etc. In general, in this work we will deal with the raw 
observed luminosities (from which we can compare our 
host galaxies to each other and to observed stellar pop¬ 
ulations free of any assumptions) at most stages of the 
analysis, converting to mass only to add physical context 
or interpret our results. 

3. RESULTS 

3.1. Near-Infrared Luminosity Distribution and Redshift 

Evolution 

In Figure [3] we plot the near-infrared luminosity (see 
above) of all galaxies in our survey as a function of red¬ 
shift. Blue curves show the magnitude that would be 
measured for several galaxy templates of various masses 

ing the resulting M 3 ,6/(i+z) (M*, 2 ) conversion function with the 
SED-fit masses and 3.6/rm magnitudes in the MODS and Ultra- 
VISTA surveys; i.e., our SED model was developed to ensure that 
we could accurately predict the average masses of MODS/Ultra- 
VISTA galaxies from just their 3.6/rm magnitudes. Indeed, single¬ 
band masses via our procedure for individual galaxies using IRAC 
photometry from the comparison field surveys show good agree¬ 
ment with masses derived from a full SED fit, with an additional 
scatter of a factor of 2—3 (equivalent to 0.3-0.5 dex). 


as a function of redshift using the templates described 
in the previous section. A few GRBs have been treated 
specially in this plot and the ensuing analysis. For the 
initial 3.6 /im observation of GRB 050401 in Cycle 4, the 
orientation of the detector resulted in a great deal of 
stray light from a bright nearby star; this target was re¬ 
observed in Cycle 8 at a better orientation but at 4.5 /im 
only; because this image is deeper (and the expected 
3.6-4.5/nn AB color close to 0.0) we use the 4.5 /im 
magnitude. The region near the afterglow position of 
GRB 090313 is strongly blended with a foreground galaxy 
and is omitted from the plot and subsequent analysis 
(pending acquiring deep ground-based data to localize 
the host centroid). GRB 080310 is affected by a diffrac¬ 
tion spike from a nearby star; we attempted to model 
and subtract it, but it is difficult to evaluate the reality 
of sources at or near the host location and a conservative 
estimate of the field depth does not reach our program 
goals, so we omit this also. 

The host galaxy luminosity distribution exhibits strong 
evolution with redshift, at least at the low -2 end: lumi¬ 
nous galaxies (e.g. those with M < —21) are common 
above 2 > 1.5 but effectively vanish from the host sam¬ 
ple below 2 < 1. To better quantify this, we calculated 
the median host magnitude as a function of redshift us¬ 
ing a moving window encompassing up to 21 objects (10 
at lower redshifts and 10 at higher redshifts). Nonde¬ 
tections are included; since the median magnitude is al¬ 
ways above the limits of typical nondetectons, our lack 
of knowledge of their actual flux does not affect the me¬ 
dian calculation. The uncertainty of this running median 
is calculated (and shown as the filled region on the dia¬ 
gram) via a simple resample-with-replacement bootstrap 
technique. 15 The resulting curves are shown in Figure [3] 
in cyan. A marked increase of 2-3 magnitudes in the me¬ 
dian luminosity is evident up to 2 ~ 1.4, at which point 
the average magnitude levels off, showing no further (sig¬ 
nificant) variation out to higher redshifts. Some of this 
behavior can be interpreted as the result of the shift of 
the effective rest-frame wavelength of the fixed filter (see 
the thin blue equal-mass curves); in particular a mod¬ 
est downturn of up to 1 mag is expected below 2 < 1 
as a result of the stellar bump moving out of the 3.6 /im 
bandpass. However, star-forming galaxy SEDs are rel¬ 
atively flat (in AB) between rest-frame wavelengths of 
0.8-2.2/rm which most of our observations correspond to, 
and the strong decline in the typical luminosity (by 2- 
3 mag) at low redshifts is too large to be interpreted as 
anything other than intrinsic evolution in the population. 
Specifically, the characteristic mass of the host popula¬ 
tion declines by approximately an order of magnitude 
from 2 ~ 1.5 to the present time. 

An alternative visualization of the redshift evolution of 
the luminosity distribution is shown in Figure U which 

15 The treatment of nondetections does affect the faint bound 
of this uncertainty envelope at z > 3.7, since we detect just above 
50% of hosts at these redshift and many replacement trials produce 
a median set by the assumed value of the limits. A maximally con¬ 
servative treatment would place the nondetections at a flux level of 
zero; in this situation the lower uncertainty on the median strictly 
becomes unbounded. For Figure [3] we randomly assign nondetec¬ 
tions at a flux level between the limiting value and zero, but we 
indicate the lower bound using a dotted line in this region to indi¬ 
cate that this is dependent on the assumption of the distribution 
of limiting fluxes. 
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plots the cumulative luminosity distribution in a series 
of redsliift bins from 0 < 2 < 0.5 to 4.5 < 2 < 6. 
Nondetections are represented arbitrarily as the limiting 
magnitude plus 1.0, and shown as a dotted line. While 
the two lowest-redshift bins show a distribution weighted 
towards low-luminosity hosts, the remaining curves are 
nearly identical. 

3.2. Dark Bursts and the Role of Extinction 

The nature of the “dark” GRB population and its 
implications for the GRB host galaxy population over¬ 
all has been a major topic in the study of GRB envi¬ 
ronments during the past decade. Systematic follow¬ 
up of the host galaxies of these events has shown that 
the hosts of the dust-obscured population show dramat¬ 
ically higher star-formation rates and stellar masses and 
redder colors than the hosts of unobscure d GRBs (e.g., 
Kriihler et al.1120111: IS vensson et al.l 120121 iPerlev et al'l 
20131 iRossi et al.l l2012t iHunt et all 120141) . with major 
implications for the nature of the GRB host population 
(and its connection to star-formation) overall. 

Given these past results, it is unsurprising that we ob¬ 
serve a strong correlation between the existence of ob¬ 
scuration affecting the afterglow and the properties of 
the host, as illustrated in Figure[5] It is noteworthy that 
GRBs occurring within the most NIR-luminous quartile 
of GRB hosts (those at M 3 6 /( 1+z ) < —22.5 mag, or a 
mass of about M* > 3 x 10 1O M 0 ) are almost exclusively 
moderately- to severely obscured. This result also em¬ 
phasizes the critical role of choosing a uniformly-selected 
sample for studies of this type: had our survey consid¬ 
ered only those events with optical afterglow-determined 
redshifts (light-colored red and green points in the main 
panel of Figure H our results would be entirely differ¬ 
ent, since nearly all of the luminous galaxies in the sam¬ 
ple hosted GRBs whose redshifts required host-galaxy 
follow-up. 

In contrast, GRBs in the least-luminous galaxies are 
only rarely obscured. Among GRBs in our sample 
at known redshift, only a single example is defini¬ 
tively dark/dusty and hosted within a galaxy with a 
luminosity-derived mass less than ~ 3 x 10 9 Mq. While 
three others at unknown redshift are in galaxies that are 
very faint and— unless they are at 2 ~ 5—likely hosted 
in similar galaxies, this is still dwarfed by the numbers 
of dusty GRBs in massive galaxies or by the numbers of 
unobscured GRBs in low-mass galaxies. A dust-obscured 
afterglow in a low-mass host is a very rare situation. 

While a detailed analysis of the coupling between 
the properties of dust seen towards the GRB after¬ 
glow ( Ay , extinction law, etc.) and host galaxy prop¬ 
erties will be reserved for future work, we note in 
passing that the nearly-ubiquitous presence of obscura¬ 
tion in the afterglows of GRBs originating from lumi¬ 
nous systems provides further evidence in support of 
the notion that dust in massive galaxies is fairly ho¬ 
mogeneous, or at least contains a large diffuse com¬ 
ponent. It also suggests that ongoing star-formation 
in dense and dusty clouds in low-mass galaxies is not 
a large contributor to cosmic star-formation and more 
broadly that low-mass galaxies harbor relatively lit¬ 
tle dust. ( See a lso our discussion of this topic in 
Perlev et al.112013L and the discussion in ISchadv et al.1 
2014 ) This is in agreement with conclusions of decreas¬ 


ing d ust content towards fainter galaxies in deep surveys 
('e.g.. iBouwens et allI2009L12014 Gastellano et al.l 120121 : 

iFinkelstein et al.1 l2012t lOesch et al.l 120131) but extends 
these results to include optically-thick, heterogeneously- 
distributed dust (which would not manifest itself in 
the galaxy colors) and to galaxies with arbitrarily low 
masses. 

We do not see obvious redshift evolution in the ten¬ 
dency for dust-obscured GRBs to inhabit more lumi¬ 
nous galaxies: all but one of the dust-obscured bursts 
at 3 < 2 < 5 lies within a host galaxy with a mass 
above the median. This suggests that even over this red¬ 
shift range, dark GRBs cannot be neglected in studies 
of host demographics. Several of these GRBs did have 
an afterglow bright enough to secure an absorption-based 
redshift measurement, however—which may suggest that 
dusty GRBs become less so between 2 ~ 2.5 and higher 
redshifts. However, since dust-obscured GRBs represent 
only a modest fraction of the total at any redshift, we do 
not yet have sufficient number statistics to address the 
issue definitely. 

3.3. Comparison to Field-Selected Star-Forming 
Galaxies 

The nature of the GRB host population at any redshift 
is affected both by the true distribution of star-formation 
within galaxies as well as any preference that the GRB 
progenitor may have for particular types of environment, 
and the primary goal of our survey is to search for and 
characterize these influences. A detailed examination of 
this topic will require completion of the multi-band sur¬ 
vey and as such will be addressed in future publications— 
however, the Spitzer observations alone are already quite 
informative, since the Spitzer magnitude serves as a good 
tracer of the stellar mass distribution and can be used 
to produce model-independent comparisons against star¬ 
forming galaxies from galaxy field surveys with Spitzer 
data. 

_Following an a pproach similar to our previous study 
(iPerlev et all20131). in Figure [3] we also plot field galaxies 
from lKaiisawa et al.l 1 2011 ). with the area of each point 
scaled proportional to the star-formation rate of each 
object. These galaxies come from a deep IF s -band se¬ 
lected catalog (MODS, the MOIRCS Deep Survey) of the 
GOODS-North field. This survey is very deep ( K s > 24.9 
AB mag at 5<r) and mass-complete to a similar level as 
is achieved by our Spitzer observations, and contains a 
large value-added catalog including dust-corrected UV 
and 24/xm star-formation rate measurements. It there¬ 
fore makes an appropriate field-survey comparison sam¬ 
ple for our study, although it is susceptible to cosmic vari¬ 
ance as a result of the small field footprint (28 arcmin 2 
for the “wide” catalog employed here). To mitigate 
against cosmic variance, we also employed the catalogs 
derived fro m deep observations of the UPS field fr om 
CANDELS (iGalametz et al.ll2013i : lSantini et al]l2015f) — 
approximately 7 times larger in area, although without 
the long-wavelength star-formation rate measurements. 
Since both these surveys have small numbers of targets 
and may be incomplete at low redshift (2 < 0.5), we 
also used the much wider (5400 arcmin 2 ) but shallower 
COSMOS/Ultra-VISTA catalogs (iMuzzin et a,l.ll2013h as 
a comparison population for the (very limited) sample of 
GRB hosts in that redshift range. 
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Fig. 5.— Same as Figure[3]but now with an apparent-magnitude cut applied to both populations in order to compare the samples only 
in the regions where both are complete. The cyan curve again shows the median magnitude of the GRB host population near that redshift 
(and lcr bootstrap uncertainty), this time only for m < 24.25 mag detections. The dashed and dotted curves show the same property 
calculated from MODS/GOODS-N and CANDELS/UDS, respectively. GRB hosts are significantly underluminous as a population at all 
redshifts but especially below z < 1.5. The red measurements show our best-fit values for M 3 . 6 /(i+ z ).thresh? a (soft) upper limit for efficient 
GRB production within a host galaxy (Table [TJ. The solid red curve shows the threshold for a metallicity of log[0/H]=8.94 converted to 
luminosity (ES- 


The contribution of undetected galaxies to the cosmic 
star-formation rate in a field survey is fundamentally 
uncertain without employing assumptions: extrapolat¬ 
ing luminosity functions or relying on theoretical mod¬ 
els. In contrast, GRBs probe star-formation regardless 
of the detectability of their host galaxies. To fairly com¬ 
pare the GRB and field galaxy luminosity distributions, 
we trim our samples at an apparent magnitude cut of 
m 3.6 < 24.25 AB, which corresponds to the complete¬ 
ness limit of both our host-galaxy survey and of MODS, 
and examine only the redshift range of 0.5 < 2 < 3, 
since even AT-band selected samples cannot be complete 
with respect to rest-frame NIR luminosity at higher red- 
shifts. (Where relevant we use the aperture-corrected 
total magnitudes, not the raw catalog magnitudes.) We 
then calculate the median absolute magnitude of the host 
galaxy distribution following the same procedure as in 
the previous section, and the SFR-weighted median ab¬ 
solute magnitude of the field galaxy sample. For this pur¬ 
pose, within MODS we use the SFRir + uv column, which 
adds UV and 24/xm star-formation rates if detected with 
Spitzer-MIPS, or otherwise uses an extinction-corrected 
UV measurement. For UDS we use the star-formation 
rate from the “14a” SED model (which employs the most 
flexible star-formation history) and for Ultra-VISTA we 
employ the SED-fit star-formation rate estimate. These 
median curves are plotted as thick solid blue (GRBs) and 


dashed black lines (MODS) in Figure [5] If GRBs were 
perfect tracers of the cosmic star-formation rate (that 
is, the GRB rate per unit star-formation was a constant 
regardless of environmental properties such as metallic¬ 
ity), then the medians of the two distributions would be 
statistically consistent with each other at every redshift. 

Large differences are observed between the actual av¬ 
erage luminosity of our GRB host sample and that of the 
(SFR-weighted) field galaxy population at low redshift: 
about 1.5 mag (a factor of four) below 2 ~ 1. The devia¬ 
tion gradually narrows with increasing redshift but does 
not completely disappear: at redshifts of 2 = 2 — 3 the 
GRB host population is still located in galaxies that are, 
on average, about 0.5 mag fainter than what would be 
predicted for a strictly uniform tracer. 

A more detailed representation of this behavior is 
shown in Figure |H1 which provides cumulative distribu¬ 
tions for the absolute magnitude of GRB hosts in vari¬ 
ous redshift ranges, compared to the absolute magnitude 
distributions weighted by star formation rate of galaxies 
in several field surveys. The MODS sample is shown 
as the black line and the UDS sample is shown as a 
gray line, except in the lowest-redshift bin where we use 
COSMOS/Ultra-VISTA. The GRB distribution is heav¬ 
ily weighted towards low-luminosity (low-mass) galaxies 
out to at least 2 ~ 1.5, and skewed more weakly towards 
higher redshifts. 
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Fig. 6. — The cumulative luminosity distribution of GRB hosts (colors) compared to the star-formation-weighted distribution of galaxies 
in GOODS-N/MODS (black) and UDS/CANDELS (gray) above the completeness limit of the survey. (For the lowest-redshift bin we 
instead compare to Ultra-VISTA.) Since field surveys are incomplete to faint galaxies the galaxy curve is anchored to the host distribution 
starting at a common limiting magnitude shown as the dotted vertical line. The curves are highly inconsistent at low redshifts with an 
almost complete dearth of GRBs in luminous and massive galaxies responsible for 30% or more of cosmic star-formation at that epoch 
(as also seen in Figure [3J. Milder inconsistency is observed at higher redshifts. No appropriate mass-limited sample exists with which to 
compare beyond z > 3.5. The dashed portion of the GRB curve incorporates nondetections in the same way as Figure[4] 


3.4. The Threshold Mass and Metallicity for GRB 
Production 

A striking feature of Figures [3] and [6] is the near-total 
absence of GRBs in massive systems at low redshift. 
Galaxies more massive than 3 x 10 10 M 0 (about half 
the mass of the Milky Way) are responsible for more 
than a quarter of the cosmic star-formation rate density 
at z < 1 but produced no GRBs in our sample, which 
is 100% complete and unbiased for galaxies in this mass 
and redshift range. In contrast, galaxies with a 3.6^im- 
derived mass only slightly less (10 9 — 1O 1O M 0 , within 
a factor of a few of the LMC) produce the majority of 
GRBs in this interval. This suggests two things. First, 
it requires that that the GRB rate is a quite sharp func¬ 
tion of environment near its critical level: e.g. the GRB 
efficiency in 10 1Q — 10 n M 0 galaxies must be roughly 
an order of magnitude lower than the efficiency within 
10 9 — 10 w Mq galaxies. However, it also indicates that 
extremely low mass is not necessarily a boon to GRB pro¬ 
duction: while some extremely faint GRB hosts do exist, 
the median mass and luminosity are still well within the 
range that is considered “normal” for star-forming galax¬ 
ies. Extremely faint galaxies are not required or even 
particularly favored for GRB prouction. 

Our result clearly rules out models in which GRBs 
are entirely uniform (unbiased) with respect to star- 
formation, as well as models in which the GRB depen¬ 
dence is relatively smooth. For ex ample, the model fa¬ 
vored by the study of lTrenti et al .1 120151 (largely on the 
basis of fitting the GRB redshift distribution, but also 
incorporating some limited host observations) suggests 
that most (~85%) GRBs at z < 1 form via a metal- 


independent channel; this would predict a large num¬ 
ber of GRBs from massive and luminous galaxies in the 
nearby universe and is in contradiction to our observa¬ 
tions. (Trenti et al. themselves also note that the lack of 
reported luminous low-z hosts is a point of tension with 
their model.) 

It is, however, in general agreement with previous 
studies of the GRB host stellar mass distrib ution at 
z < 1. For example, Figure 3 of iBoissier et al.l 120131 ) 
also shows an order-of-magnitude drop in GRB efficiency 
between a host stella r mass of 10 9 6 and 10 10 4 Mq, while 
IVergani et al.l (120151 ) find no GRB hosts at > 1O 1O M 0 in 
their complete z < 1 sample (requiring strong supression 
in this range) but measure a variation in efficiency of only 
a factor of a few or less below 1O 95 M 0 (their Figure 4). 
It is also qualit atively consistent with the conclusions of 
ISvensson et alJ ( 2010 ). who compared GRB and cc-SN 
host properties. 

The GRB dependence on host environment is widely 
speculated to intrinsically be a metallicity dependence, 
in which metal-poor stars form GRBs readily but metal- 
rich stars much less or not at all. Due to the com¬ 
plex correlations between different properties in galax¬ 
ies, it is difficult to be certain from photometric data 
alone whether or not metallicity (and not some other 
parameter that also correlates with mass/luminosity) is 
actually responsible for the mass-dependent GRB rate 
at low redshift—other possibilities invoke variations in 
the IMF/binarity associated with starburst intensity, or 
stellar mergers associated with dynamical interactions 
found only Jn very d ense/intense star-forming regions 
(jvan den Heuvel fe Portegies Zwartll2013l) . Direct spec- 
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troscopic metallicities of all of the massive galaxies in our 
sample would provide a means to evaluate this directly, 
but we do not (yet) have these observations for our full 
sample—and even with spectroscopy in hand, ambigui¬ 
ties associated with inhomogeneous conditions in these 
galaxies can introduce significant complications in inter¬ 
pretation. 

However, our sample offers an opportunity to investi¬ 
gate the metallicity-limit hypothesis in a different way. 
The mass- metallicity relation evolves strongly with red- 
shift (e. g.. lErb et al.l l2006lh leading to a natural pre¬ 
diction (iKocevski et al.lT2009D that the threshold stellar 
mass (and NIR luminosity) should rise with redshift, and 
do so in a specific way. Previous studies (iKruhler et al.1 
1201 It iPerlev et al.11201311 have hinted that this effect is 
present but have been unable to make quantitative tests 
of it due to their incomplete and potentially biased red- 
shift sampling for z > 1 objects. With our large, com¬ 
plete sample we are in a position to quantitatively exam¬ 
ine this behavior for the first time. 

With a simple metallicity-cut model in mind, we use 
our redshift-subdivided luminosity distributions (Figure 
© to test whether the increasing stellar mass with red- 
shift that is apparent in our sample is indeed consistent 
being the product of a redshift-dependent mass metallic- 
ity relation and a redshift-invariant metallicity threshold 
for GRB production. For each redshift interval, we re¬ 
sample the field galaxy magnitude distribution such that 
it is uniform with respect to star-formation (binning in¬ 
dividual galaxies of similar magnitudes to produce a new 
distribution where each magnitude increment has con¬ 
stant total SFR, mirroring the way GRBs are expected 
to select galaxies) and perform a two-sided K-S test be¬ 
tween this and the GRB host luminosity distribution. 
This procedure is then repeated for a large number of 
possible maximum-luminosity thresholds spanning the 
range of the host sample, recalculating the K-S p-value 
each time. We then estimate the threshold level, and 
its upper/lower uncertainties, by marginalization. These 
measurements and uncertainties are presented in Table 

m 

Despite the large size of our survey, we remain strongly 
limited by the number counts of the GRB host popu¬ 
lation in each bin and the uncertainty on the thresh¬ 
old measurement is large (typically 0.5-1.0 mag). How¬ 
ever, at all redshift bins at which we can carry out this 
analysis, a model with no luminosity cut is ruled out 
(95 — 99.99% confidence, depending on redshift). We do, 
however, find a good match to the distributions at every 
redshift if a luminosity cut is applied to the field-galaxy 
sample. 

Our measured luminosity thresholds (and uncertain¬ 
ties) can then be translated into stellar masses using the 
same procedure described in 112.41 From there, they can 
be converted into metallicity values using the analytic 
expressions for th e red shift-dependent mass-metallicity 
relation of lZahid et, al.l (|2014l ). 

If the GRB efficiency is controlled only by metallicity 
(and not other factors that are varying over the Uni¬ 
verse’s history, such as the average specific SFR) we 
would expect to measure a consistent metallicity thresh¬ 
old at each redshift. Indeed, we find that all of our ob¬ 
servations are consistent with redshift-independent GRB 
metallicity threshold of 12 + log[0/H] = 8.94 on the 


KK04 system with a scatter of only a = 0.04 among 
the diff erent redshift bins. Co nverting using the rela¬ 
tions of lKewlev fe Ellisonll2008l this is equivale nt to 12 
+ log[Q/H] = 8.64 on the lPettini fe Paeell (120041) (PP04) 
system. In eith er case th is is close to the Solar value 
(~8.7. lAsnlund et, al.ll2009h . 

Given that the original Zahid et al. relation was de¬ 
rived from a fit using only z < 1.7 data the extension 
of their analytic prescription to higher redshifts is purely 
an extrapolation, and it is somewhat surprising that we 
continue to see agreement between our measured thresh¬ 
old and their mass-metallicity-redshift relation even at 
2 = 2.25 and (to some extent) at 2 = 3. The actual mass- 
metallicity relation at higher redshi fts has been i nvesti¬ 
gated directly by other surveys (e.g..lSavaglio et al.l[200.5l 
lErb et all 120061: iMaiolino et al.ll2008l) ls butT remains un¬ 
certain, in part because the ionization properties of high- 
rcdshift galaxies appea r to be quite different than those 
at low er redshifts ('e.g.. lSteidel et al.1l20l4 ISanders et al.l 
120151 ), putting the calibration of popular diagnostics and 
prescriptions developed at z ~ 0 into doubt at these 
epochs. Our observations may actually help shed some 
light on this topic: given the good consistency between 
the GRB host luminosity function with the evolving 
mass-metallicity relation from z ~ 0.5 to z ~ 2.0 with a 
fixed, approximately Solar metallicity limit, the contin¬ 
ued avoidance of GRBs from massive (~ 10 11 M e ) galax¬ 
ies out to z ~ 3 suggests that massive systems are metal- 
enriched at or above at least the Solar value (and would 
favor the use of abundance diagnostics/calibrations con¬ 
sistent with this evaluation). If so, this would be an in¬ 
teresting reversal: while the metallicity sensitivity of the 
GRB progenitor has long been seen as an obstacle for its 
use to measure the cosmic star-formation rate density, 
this property may in the end turn out to be valuable by 
making GRBs a independent tracer of chemical evolu¬ 
tion. 

Turning to low redshifts, although our inferred cut 
value provides a good explanation for the luminosity dis- 
trubution of galaxies over the full range of redshifts well- 
probed by our survey, it is somewhat higher than the 
effective metallicity cuts measured b y spectrosco pi c sur- 
veys at low redshift. For example, iModiaz et al.l 120081 
infer a cutoff value of ^8.5 (on the iKewlev fc Donital 
120021 abundance scale, equivalent to 8.67 on KK04) based 
o n fou r objects at z ~ 0.1; the updated analysis of 
I Graham fc Fruchteil (2013! finds reasonable agreement 
with this cut value for the bulk of their population (al¬ 
though fully 5 of their 14 hosts have measured m etal- 
licities above this.) IWolf fc Podsiadlowskii (120071 ) find 
a cutoff value of 8.7 ± 0.3 (KK04) using a photometric 
analysis of published low -z hosts. Given the large uncer¬ 
tainties on these earlier values (associated with the small 
sample size) they might not be inconsistent with our 
new measurement. In addition, it must be remembered 
that the sample of local, low-luminosity GRBs has prop¬ 
erties quite different from the cosmological population 
observed at high redshift (e.g., >5 orders of magnitude 
in Eiso and in observed volumetric rate; low-luminosity 
GRBs are also quasi-spherical and probably not jetted; 

16 It has also been studied using GRB hosts themselves, compar- 
in g line-of-sight ab sorption metallicities to Spitzer stellar masses 
(ILaskar et a.I.lLOl III . 















































12 


Perley et al. 


TABLE 1 

Redshift-Dependent Cutoff Values 



z a 

M 3.6/(l+z ) b 


12+log[0/H] d 

Al3.6/(l+z),8.94 e 

/SF f 

Pbest^ 

P8.94 h 

Poo 1 

0 .1- 

-0.5 




-20.30 





0.5- 

-1.0 

—21 33 —0,31 

^x.oo_|_Q 34 

9-82+0J5 

8-91 !g;8S 

-21.53 

0.33 

0.71 

0.43 

<0.0001 

1 .0- 

-1.5 

—23 17 —1,23 
zo - ± ‘+0.97 

10 70^" 0-57 

1U. /O_o 52 

9.05!°!? 

-22.23 

0.43 

0.24 

0.06 

0.044 

1.5- 

-2.0 

—22 41 ~ 0-36 

ZZ.H:±_|_ 0 .54 

10.341°;*° 

8 Qi +0-05 

-22.57 

0.50 

0.13 

0.12 

0.003 

2 .0- 

-2.5 

r\n oi — 0.38 

zo.oi_|_Q 62 

10.831°;*° 

8 99+ 0 - 04 

o.yy _ 0 09 

-22.92 

0.62 

0.70 

0.19 

0.017 

2.5- 

-3.5 

—22 97 “ 0-27 
' +0.48 

10-64!“;“ 

o OO+0.04 

o-oo_0.08 

-23.39 

0.77 

0.87 

0.16 

0.012 


a Redshift range 

b Inferred AB absolute magnitude threshold and 15-85% confidence uncertainty interval, measured from the 
data. Confidence interval does not include uncertainty due to cosmic variance in the field survey. 
c Stellar mass corresponding to column (b). 

d Oxygen abund ance corresponding to column (c), converted using the analytic mass-metallicity relation of 
IZahid et al.1 (120141 1. 

e Absolute magnitude threshold predicted, assuming a metallicity cutoff of 12+log[0/H] = 8.94 at all redshifts. 
f Fraction of cosmic star-formation at that epoch which readily produces GRBs. 

s p value calculated from a two-sided K-S test betwen the GRB host magnitude distribution and the SFR- 
weighted galaxy distribution clipped at the best-fit magnitude threshold. 
h p value if clipped at the magnitude threshold corresponding to 12+log[0/H] = 8.94. 

1 p value if not clipped at all (if GRBs are assumed to be “unbiased” tracers of SFR). 


iBromberg et al.l 1201 If) . Most low-redshift samples are 
also still based primarily on pre-Swift GRBs whose selec¬ 
tion is particularly heterogeneous and difficult to model; 
future studies of low -0 events from the more uniform, 
better-understood Swift population should help better 
constrain the true nature of the low-redshift population. 

3.5. Peering Below the Surface: The Fraction of Star 
Formation in IRAC-Undetected Galaxies 

If the GRB efficiency is close to uniform in galaxies 
with luminosities/metallicities below our inferred thresh¬ 
old, our observations can be used to provide an inde¬ 
pendent means of estimating what fraction of all star- 
formation occurs in the least luminous and lowest-mass 
galaxies at each redshift. At the completeness level of 
the MODS survey (approximately 7713.6 = 24.25, slightly 
shallower than the level achieved our host galaxy sur¬ 
vey) the star-forming fraction in fainter galaxies can be 
read directly from our plots in Figure [7] as the fraction 
where the completeness magnitude intersects the GRB 
host curve. Given the depth of the survey, this frac¬ 
tion is unsurprisingly quite low at 0 < 2: 10-20% in 
every well-sampled redshift bin. At higher redshifts the 
implied survey-detectability fraction drops significantly: 
the fraction of faint hosts is ~30% at 0 = 2.25 and ~50% 
at 0 = 3. Our survey goes deeper than MODS at 0 ~ 3 
(77iii m ~ 25) and gives an even tighter constraint to this 
depth; only about 30% of star-formation is in galaxies 
fainter than that level (about 3 x10 9 Mq). 

Extending this constraint to even higher redshift (0 > 
3.5) is contingent on the assumption that the GRB host 
distribution probes the entire galaxy population (i.e., 
that even the most massive galaxies have metallicities 
below our inferred threshold), an assessment we cannot 
test directly because no mass-complete field-galaxy sam¬ 
ples exist at these redshifts. 1 ' However, at the depths 

17 Two other manuscripts were recently submitted on this topic, 
using the rest-frame UV lu minosit y distribution. They come to 
different conclusions: ISchulze et al.ll2015l show evidence for a GRB 


achieved on our deep high -0 fields ( 7773.6 ~ 25.5 mag, 
corresponding to ~few xlO 9 M 0 ) we continue to de¬ 
tect approximately half the GRB population out to the 
redshift limit of the sample (0 ~ 5.5), indicating that 
about half of star-formation occurs in galaxies less lumi¬ 
nous/massive than this threshold. If GRBs continue to 
avoid high-mass galaxies even in this redshift range, the 
implied fraction of undetectable star-formation would be 
less than this. 

3.6. Correcting the GRB-Inferred Cosmic 
Star-Formation Rate History 

Figure [7] also allows us to measure directly the frac¬ 
tion of star formation at each epoch that does vs. does 
not produce GRBs (except possibly at a much reduced 
rate)—this is simply the fraction of the plot below vs. 
above (respectively) the horizontal dotted line. In the 
last column of Table [T| we provide this measurement, 
which is ~ 0.85 at 0 ~ 3 (indicating that GRBs track 
star-formation well in all galaxies except for those re¬ 
sponsible for 15% of cosmic star-formation) but drops 
to 0.5 at 0 ~ 1.5 (indicating that half or more of the 
star-formation in the Universe does not produce GRBs, 
except perhaps in unusual circumstances) and even lower 
at lower redshifts. 

This allows us to correct the cosmic GRB rate esti¬ 
mated in Paper I to examine whether or not the same 
model that explains the evolution of the luminosity dis¬ 
tribution can also explain the redshift distribution and 
in particular the apparent excess of GRBs at high -0 (or 
equivalently, the deficiency at low- 0 ). We use a polyno¬ 
mial fit to interpolate between the redshift bins for the 
values in Table Q] and simply divide each bin by the re¬ 
sulting value of /sf- 

The result is plotted in Figure [8] (left panels). The cor¬ 
rection procedure greatly reduces the discrepancy; while 

host distribution that may be skewed towards fainter hosts at 
z ~ 4 (but not at z ~ 2 — 3) compared to a uniform tracer, while 
I Greiner et al .1 [20151 argue that the GRB host luminosity distribu¬ 
tion follows the expected distribution. 
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Fig. 7. — Similar to Figure[6] but u nder a model in which GRB production is stifled above a certain luminosity threshold. The threshold 
is determined via our procedure in by marginalization using a K-S analysis, smoothed in redshift via a polynomial fit: this value 
is shown as the right vertical dotted line in each panel. The threshold is also translated to the fraction of star-formation traced well by 
GRBs, indicated by the horizontal dashed line. GRBs are poor tracers of star-formation at low redshifts (nearly all GRBs at 2 < 1.5 were 
produced by galaxies responsible for 50% or less of cosmic star formation at that epoch). At higher redshifts they are much better tracers, 
but are still not “perfect” (the ~15% of star-formation in the most luminous galaxies still produce few to no GRBs). 
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Fig. 8 .— The cosmic rate density of GRBs compared to that of star-formation. The top panels show the raw, uncorrected GRB rate: 
the plot on the left is based on the SHOALS sample and uses the same data and calculation procedure as Paper I (based on S and 
£ is o); the plot on the right uses all Swift GRBs from 2005-2014 (a larger but poten tially biased sample) and uses the more traditional 
/av and Li so ,av Error bars show the 10-90 % binomial confidence interval (flC amero nl 1201 111. An “excess” in the GRB rate relative to the 
galaxy-inferred star-formation rate density (Madau & Dickinson 2014; Robertson et al. 2015) is seen at high- 2 . In the bottom panels we 
have corrected these rates based on the observed tendency for GRBs to avoid luminous galaxies at low- 2 . Applying this correction produces 
good consistency in the respective rates, with the exception of the 2 = 8 — 10 flux-limited bin. 
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the z > 2 data points are still above the SFRD curves 
the excess is less than a factor of two and always less 
than 2cr significance. 

Since our GRB rate derivation procedure (using E[ so 
from a fluence-limited sample) differs from previous 
work, in the right panels of Figure[8]we also show the rate 
derived in the more conventional fashion using the aver¬ 
age flux/luminosity of each burst, F = S/T 90 (we con¬ 
tinue use a fc-correction to the 45-450 keV band instead 
of bolometric corrections). Results are nearly identical, 
with the exception that this procedure recovers two z > 8 
bursts that were dropped by our f l uence cut (090423 
and 090429B; [S alvatcr raet al.ll2009t iTanvir et all 120091 : 
ICucchiara et al.l 12011 1. which would imply a very high 
rate at this epoch. This may imply different behavior 
entering the reionization era (in the burst rate or per¬ 
haps in burst properties—both events were remarkably 
short after taking into account time dilation) but may 
also be a fluke. More searches for very high -2 GRBs will 
be needed to investigate this. 

Otherwise, our analysis provides further support to the 
notion that the entire GRB rate behavior—both as a 
function of galaxy mass and as a function of redshift 
can be explained in terms of the chemical history of the 
universe. Above z > 3 GRBs track star-formation; be¬ 
low this point massive galaxies build up metals and cease 
GRB production, shifting the host luminosity distribu¬ 
tion to fainter galaxies and reducing the cosmic GRB 
rate. 

4. DISCUSSION AND CONCLUSIONS 

In this paper we have demonstrated the signatures of 
cosmic evolution in the physical properties of the GRB 
host population across a vast redshift range from z ~ 0.3 
out to z ~ 6, only 1 Gyr from the Big Bang. The me¬ 
dian host NIR luminosity does not evolve much between 
z ~ 5 and z ~ 1.5, but at lower redshifts (z < 1.5) the 
average luminosity drops by over a factor of 10. The 
2 < 1.5 host population is dramatically fainter than 
would be expected for an population that traces star- 
formation uniformly, and our observations require that 
the GRB rate per unit star-formation is greatly reduced 
in massive galaxies compared to low-mass galaxies. At 
higher redshifts of 1.5 < 2 < 3.5, the population is 
also slightly underluminous compared to a fully uniform 
tracer, but the difference is much more modest. Consis¬ 
tent conclusions have been reached in recent years based 
on optical (rest-fram e UV) data fe.g.. ITanvir et al.ll2012t 
iGreiner et alJ l2015h : our study links these efforts using 
a single meth odo l ogy across nearly all of cosmic history 
(see also lSchulze et al.ll2015 ). 

We emphasize that our target population was selected 
in an unbiased manner and our redshift completeness is 
very high (above 90%), so our conclusions involve no sys¬ 
tematic error associated with sample selection and min¬ 
imal systematic error associated with redshift measure¬ 
ment. In principle, since the galaxies that have escaped 
redshift measurement are all quite faint, or since we could 
in principle have misidentified a bright foreground sys¬ 
tem with the host galaxy in a small number of cases, the 
host galaxy population might actually be slightly fainter 
than we infer. But even this would leave our statistical 
comparion against star-forming galaxies effectively unal¬ 
tered: the unknown -2 hosts are fainter than the MODS 


galaxies against which our comparison is based (and so 
would not enter into the comparison); likewise, any hypo¬ 
thetical faint background hosts (which would be rare in 
any case, not re presenting more tha n a few percent of the 
sample; see e.g. iBloom et al.ll2002l or iPerlev et alJ l2013) 
would have to be undetected in our Spitzer/optical/HST 
imaging and would likewise be dropped from the com¬ 
parison. In any case, there is no possibility that the 
host galaxy population is brighter than we have reported, 
which firmly rules out any model for which the GRB rate 
is independent of environment. 

Our results illustrate the importance of an unbiased 
approach similar to our own when addressing GRB de¬ 
mographics at high redshift. Had dark bursts been omit¬ 
ted from our survey (or no campaign had been under¬ 
taken to measure their redshifts) nearly all objects in 
galaxies with a mass of 10 1o Mq or above would have 
been dropped, and the mass/metallicity threshold we 
would have inferred would have been correspondingly 
much lower than the near-Solar value we estimate from 
this study. Or, had we drawn our sample from the gen¬ 
eral Swift pool without imposing the optimized observ¬ 
ability cuts discussed in Paper I to produce a high ini¬ 
tial redshift completeness, a large fraction of the GRBs 
in fainter hosts would not be able to have their redshifts 
measured, producing an apparent host population biased 
towards brighter galaxies. These types of biases have be¬ 
deviled previous work on this topic, leading to the con¬ 
fusion regarding the degree to which GRBs do or do not 
track the star-for mation rate. (See, e.g.. lSavaglioll2012l or 
iHunt et al.ll2014j for recent, contrary views of the GRB 
host mass distribution using samples that were not drawn 
in an unbiased manner.) We hope that our unbiased se¬ 
lection, large sample, and uniform observational strategy 
will clear this point up unambiguously and pave the way 
towards effectively using GRBs to probe star-formation 
and galaxy evolution within the (large) regime in which 
they do seem to track predictions well, and perhaps even 
using their metallicity dependence as a tool for under¬ 
standing chemical evolution. 

We hypothesized that the behavior seen in our 
sample—a strong shift with decreasing redshift towards 
a lower-mass host population that deviates from the gen¬ 
eral population of star-forming galaxies more strongly— 
was the result of time-evolution in the mass-metallicity 
relation in a ma nner similar to what wa s previously pre¬ 
dicted by, e.g., iKocevski et alJ ( 2009 ). Using a sim¬ 
ple model in which the GRB efficiency is constant at 
low metallicity but falls sharply above a maximum- 
metallicity threshold, we found good agreement for 
a redshift-independent threshold metallicity of 12 + 
log[0/H] = 8.94 (on the KK04 abundance scale, equiva¬ 
lent to 8.64 on the PP04 scale)-close to the Solar value. 
The same model provides a good explanation for the 
GRB rate history, largely eliminating the previously- 
reported GRB rate “excess” at high -2 and producing a 
consistent picture (within a factor of 2) of the cosmic 
star-formation history with the results of high-redshift 
field surveys. A hint of a possible upturn at 2 ~ 8 will 
require future investigation. 

The metallicity threshold we infer is somewhat, 
but not dramatically, higher than what has been re¬ 
ported by prev ious GRB host surveys (and assumed 
by IKocevski et al.l 12009 ). On the other hand, it is 
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unambiguously much higher than the metallicity up¬ 
per li mits predicted by single-sta r theoretical models 
(e.g.. iMacFadven fe Wooslevl Il999t iHirschi et al.l 120051 : 
lYoon fc Langeil l2005t lLanger fc Normanl 120061) — likely 
indicating th at the progen it or is a binary system of some 
sort (see also lTrenti et alJl2015 ). Future models for the 
progenitor formation process will need to consider why 
the GRB production rate appears to exhibit a sharp 
metallicity threshold but does so at a value much higher 
than originally predicted. 

We emphasize that the suppression of the GRB rate 
in galaxies above our luminosity cut need not be to¬ 
tal: GR Bs in metal-rich environments are known to ex¬ 
ist (e.g.jLevesaue et al.ll2010HGraham fc Fruchteil[2013l 
iKriihler et al.1120151 : lHashimoto et al.l 1201.51) . and even if 
the metal avoidance was total we would expect chemical 
inhomogeneity to smear the threshold somewhat (jNiino 
Eon . Indeed, we observe a few hosts with luminosities 
above our best-fit threshold at each redsliift range—but 
always fairly close to the threshold itself, and in limited 
numbers that suggest that the rate in this type of host 
must be lower by about an order of magnitude comped 
to hosts below the threshold. 

If our interpretation that the GRB efficiency is nearly 
constant below this threshold is correct, our results can 
also be used to constrain the importance of low-mass 
galaxies to cosmic star-formation. We detect a large 
majority of the targeted hosts at low to moderate red- 
shift (90% at z ~ 1 and 70% at z ~ 2 — 3), consis¬ 
tent with the notion that most star-formation occurs in 
galaxies in the range of masses detectable to Spitzer at 
these distances (approximately 2 > 1O 9 M 0 )—not in ex¬ 
tremely low-mass dwarfs or outside galaxies, consistent 
with the present understanding of star and galaxy for¬ 
mation (but in contradiction to a recent stu dy of the 
unresolved optical background: IZemcov et al.ll2014l ). We 
continue to detect about half of the targets in the sam¬ 
ple even out to z ~ 5.5, indicating that we have not yet 
reached the epoch of cosmic history where most of cos¬ 
mic star-formation shifts to small galaxies that are still 
being formed. Deep rest-frame UV observations of the 
(not-unbiased) sample of known GRB hosts at 2 > 5.5 
have not yet detected any targets down to very faint lev¬ 
els (jTanvir et al.ll2012[ ), which may suggest that the tran¬ 
sition to low-mass galaxies occurs close to this epoch. 

While our observations support the leading paradigm 
in which metallicity is the primary factor controlling 
the rate of GRB production as a fraction of star for¬ 
mation, they do not yet establish that it is the sole 
driver. The consistency of the host mass population 
with a metallicity-threshold model at many different 
redshifts and the similar consistency with the cosmic 
star-formation rate both suggest that other factors can¬ 
not be predominant unless they correlate with metallic¬ 
ity in very specific ways to mimic a metallicity effect. 
However, other influences remain possible. We (e.g., 
iPerlev et al.1 120131 l2015al f have previously noted that 
even the intermediate-mass host population at 2 ~ 1 — 2 
shows other peculiarities compared to what are thought 
to be typical star-forming galaxies: their colors are bluer, 
their specific star-formation rates are higher, and they 
are m ore likely to be detected at radio wavelengths. 
iKellv et al.l (1201411 has noted that GRB hosts at z ~ 0.5 
appear to be unusually compact compared to supernova 


host galaxies of similar mass. Finally, the cutoff value we 
infer from our sample at z = 0.5 —2.0 is slightly in excess 
of what has previously been inferred from lower-z hosts; 
this could simply reflect the limited size and selection 
of low-redshift samples but could also indicate the influ¬ 
ence of a second parameter (sSFR, IMF, etc.) We are 
currently collecting a large library of multi-band optical, 
NIR, and radio observations to examine these possibili¬ 
ties in detail and will report the results in future papers 
associated with the survey. 
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TABLE 2 

SPITZER 3.6/iM HOST GALAXY PHOTOMETRY 


GRB 

z 

RA a 

dec a 

m 3.6 b 

M 3 .6/(1+z) C 

logio M, d 

lexp 6 

Prog ID f 

050128 

<5.5 

14:38:17.71 

-34:45:54.6 

>25.00 



5400 

90062 

050315 

1.9500 

20:25:54.16 

-42:36:02.1 

23.43+0.08 

-21.29 

9.77 

5400 

90062 

050318 

1.4436 

03:18:51.01 

-46:23:44.1 

>25.22 

> -18.90 

< 8.63 

5400 

90062 

050319 

3.2425 

10:16:47.91 

+43:32:54.1 

24.64+0.12 

-21.02 

9.69 

7200 

40599 

050401 

2.8983 

16:31:28.79 

+02:11:14.0 

24.29+0.39 

-20.89 

9.61 

3700 

40599 

050525A 

0.606 

18:32:32.58 

+26:20:22.4 

>24.61 

> -17.64 

< 8.11 

900 

3653 

050726 

<3.5 

13:20:11.95 

-32:03:51.2 

>24.44 



5400 

90062 

050730 

3.9693 

14:08:17.10 

-03:46:17.6 

>25.47 

> -20.54 

< 9.46 

7400 

40599 

050802 

1.7102 

14:37:05.84 

+27:47:12.3 

24.86+0.29 

-19.60 

9.00 

5400 

90062 

050803 

4 q+0.6 

23:22:37.85 

+05:47:08.8 

25.49+0.42 

-20.66 

9.53 

5400 

90062 

050814 

5.3 

17:36:45.39 

+46:20:21.6 

24.48+0.11 

-22.02 

10.22 

6900 

40599 

050820A 

2.6147 

22:29:38.11 

+ 19:33:36.7 

24.85+0.28 

-20.42 

9.38 

3600 

40599 

050822 

1.434 

03:24:27.21 

-46:01:59.7 

22.60+0.08 

-21.50 

9.85 

540 

272 

050904 

6.295 

00:54:50.88 

+ 14:05:09.5 

>24.98 

> -21.79 

< 10.07 

7200 

20000 

050922B 

4 Q+ 0 - 3 

00:23:13.30 

-05:36:17.8 

24.59+0.22 

-21.78 

10.08 

5400 

90062 

050922C 

2.1995 

21:09:33.08 

-08:45:30.2 

>25.34 

> -19.61 

< 9.01 

3700 

40599 

051001 

2.4296 

23:23:48.77 

-31:31:24.0 

23.20+0.08 

-21.94 

10.08 

5400 

90062 

051006 

1.059 

07:23:14.09 

+09:30:20.1 

20.50+0.01 

-22.97 

10.58 

2700 

90062 

060115 

3.5328 

03:36:08.33 

+ 17:20:42.6 

25.40+0.31 

-20.41 

9.43 

7400 

40599 

060202 

0.785 

02:23:22.91 

+38:23:04.3 

22.59+0.04 

-20.23 

9.32 

5400 

90062 

060204B 

2.3393 

14:07:14.93 

+27:40:36.3 

22.74+0.05 

-22.32 

10.29 

5400 

90062 

060210 

3.9122 

03:50:57.37 

+27:01:34.0 

23.37+0.06 

-22.62 

10.46 

7200 

40599 

060218 

0.0331 

03:21:39.69 

+ 16:52:01.8 

20.77+0.04 

-15.01 

7.20 

700 

90062 

060306 

1.559 

02:44:22.91 

-02:08:54.2 

21.54+0.02 

-22.73 

10.50 

5400 

90062 

060502A 

1.5026 

16:03:42.63 

+66:36:02.6 

>24.67 

> -19.53 

< 8.96 

5400 

90062 

060510B 

4.9 

15:56:29.48 

+78:34:12.1 

25.15+0.17 

-21.22 

9.81 

13000 

20000 

060522 

5.11 

21:31:44.86 

+02:53:09.7 

>26.37 

> -20.07 

< 9.31 

13800 

20000 

060526 

3.2213 

15:31:18.34 

+00:17:04.9 

25.48+0.29 

-20.17 

9.30 

7200 

40599 

060607A 

3.0749 

21:58:50.4 

-22:29:47.2 

>25.05 

> -20.52 

< 9.45 

7200 

40599 

060707 

3.4240 

23:48:19.06 

-17:54:17.7 

24.10+0.13 

-21.66 

9.99 

6180 

40599 

060714 

2.7108 

15:11:26.41 

-06:33:58.2 

25.23+0.26 

-20.11 

9.25 

3600 

40599 

060719 

1.5320 

01:13:43.71 

-48:22:51.2 

22.77+0.06 

-21.47 

9.84 

1600 

70036 

060729 

0.5428 

06:21:31.78 

-62:22:12.0 

24.06+0.19 

-17.95 

8.31 

1600 

90062 

060814 

1.9229 

14:45:21.32 

+20:35:10.5 

21.44+0.03 

-23.25 

10.82 

900 

70036 

060908 

1.8836 

02:07:18.41 

+00:20:31.3 

24.73+0.20 

-19.92 

9.15 

3600 

40599 

060912A 

0.937 

00:21:08.12 

+20:58:17.7 

21.56+0.03 

-21.65 

9.91 

2700 

90062 

060927 

5.467 

21:58:12.01 

+05:21:48.6 

>25.77 

> -20.78 

< 9.63 

7500 

40599 

061007 

1.2622 

03:05:19.59 

-50:30:02.3 

23.71+0.15 

-20.13 

9.22 

5400 

90062 

061021 

0.3463 

09:40:36.13 

-21:57:04.9 

>24.78 

> -16.22 

< 7.63 

1600 

90062 

061110A 

0.7578 

22:25:09.84 

-02:15:31.4 

24.37+0.29 

-18.38 

8.41 

1600 

90062 

061110B 

3.4344 

21:35:40.40 

+06:52:33.9 

>25.24 

> -20.52 

< 9.47 

7400 

40599 

061121 

1.3145 

09:48:54.54 

-13:11:43.2 

21.50+0.01 

-22.43 

10.31 

5400 

90062 

061202 

2.253 

07:02:05.68 

-74:41:55.3 

23.24+0.06 

-21.75 

9.99 

5400 

90062 

061222A 

2.088 

23:53:03.40 

+46:31:58.8 

24.06+0.08 

-20.79 

9.55 

18500 

30000 

070110 

2.3521 

00:03:39.24 

-52:58:27.3 

>25.13 

> -19.95 

< 9.16 

3600 

40599 

070129 

2.3384 

02:28:00.91 

+ 11:41:04.2 

23.01+0.03 

-22.05 

10.15 

6200 

40598 

070223 

1.6295 

10:13:48.40 

+43:08:00.6 

23.59+0.08 

-20.77 

9.53 

5400 

90062 

070306 

1.4959 

09:52:23.31 

+ 10:28:55.4 

21.40+0.03 

-22.79 

10.53 

900 

70036 

070318 

0.840 

03:13:56.77 

-42:56:46.2 

>24.30 

> -18.67 

< 8.54 

1600 

80153 

070328 

2.0627 

04:20:27.65 

-34:04:00.5 

23.12+0.06 

-21.71 

9.97 

5400 

90062 

070419B 

1.9588 

21:02:49.78 

-31:15:49.1 

23.30+0.09 

-21.43 

9.84 

5400 

90062 

070508 

0.82 

20:51:11.72 

-78:23:04.8 

22.46+0.07 

-20.46 

9.41 

1600 

70036 

070521 

2.0865 

16:10:38.68 

+30:15:22.8 

21.93+0.04 

-22.92 

10.63 

1600 

70036 

070621 

<5.5 

21:35:10.05 

-24:49:02.5 

>24.97 



5400 

90062 

070721B 

3.6298 

02:12:32.96 

-02:11:40.8 

>25.47 

> -20.39 

< 9.42 

10800 

80054 

070808 

~1.35±0.85 

00:27:03.30 

+01:10:35.4 

23.57+0.10 

-20.41 

9.35 

5400 

90062 

071020 

2.1462 

07:58:39.78 

+32:51:39.6 

>24.36 

> -20.54 

< 9.43 

3600 

80054 

071021 

2.4520 

22:42:34.32 

+23:43:06.2 

22.10+0.05 

-23.05 

10.68 

1600 

70036 

071025 

00 
+ 
o o 

^ 4^ 

23:40:17.07 

+31:46:42.6 

24.32+0.15 

-22.01 

10.18 

4500 

80153 

071112C 

0.8227 

02:36:50.97 

+28:22:16.5 

23.59+0.11 

-19.34 

8.89 

2700 

90062 

080205 

9 79+O.24 
‘ z —0.74 

06:33:00.69 

+62:47:32.0 

23.79+0.13 

-21.55 

9.91 

5400 

90062 

080207 

2.0858 

13:50:02.97 

+07:30:07.3 

20.98+0.02 

-23.87 

11.11 

1600 

50562 

080210 

2.6419 

16:45:04.01 

+ 13:49:35.6 

>24.62 

> -20.67 

< 9.50 

1600 

80153 

080310 

2.4274 

14:40:13.80 

-00:10:30.7 

23.84+0.25 

-21.29 

9.78 

3500 

80054 

080319A 

2.0265 

13:45:20.01 

+44:04:48.4 

22.63+0.03 

-22.16 

10.21 

5400 

90062 

080319B 

0.9382 

14:31:41.00 

+36:18:08.6 

24.57+0.22 

-18.64 

8.50 

3600 

11116 

080325 

1.78 

18:31:34.22 

+36:31:23.7 

21.09+0.02 

-23.45 

10.91 

1600 

70036 

080411 

1.0301 

02:31:55.19 

-71:18:07.3 

>24.66 

> -18.75 

< 8.54 

2700 

90062 

080413A 

2.4330 

19:09:11.75 

-27:40:40.4 

>24.15 

> -20.99 

< 9.64 

3600 

80054 

080413B 

1.1014 

21:44:34.67 

-19:58:52.3 

23.22+0.09 

-20.33 

9.30 

2700 

90062 

080430 

0.767 

11:01:14.71 

+51:41:08.4 

24.55+0.34 

-18.22 

8.33 

1600 

90062 

080603B 

2.6892 

11:46:07.67 

+68:03:39.8 

25.43+0.24 

-19.89 

9.15 

10800 

11116 

080605 

1.6403 

17:28:30.04 

+04:00:55.7 

21.61+0.10 

-22.77 

10.53 

1600 

80153 

080607 

3.0368 

12:59:47.24 

+ 15:55:10.8 

22.96+0.06 

-22.58 

10.45 

4500 

70036 



SHOALS II: NIR Luminosity Distribution 


TABLE 2 

SPITZER 3.6/iM HOST GALAXY PHOTOMETRY 


080710 

0.8454 

00:33:05.63 

+19:30:05.4 

>23.52 

> -19.47 

< 8.95 

2700 

90062 

080721 

2.5914 

14:57:55.83 

-11:43:24.6 

>24.30 

> -20.96 

< 9.63 

3600 

80054 

080804 

2.2045 

21:54:40.18 

-53:11:04.8 

24.76+0.23 

-20.19 

9.28 

3600 

80054 

080805 

1.5042 

20:56:53.44 

-62:26:39.4 

22.70+0.08 

-21.50 

9.86 

1600 

80153 

080810 

3.3604 

23:47:10.49 

+00:19:11.5 

23.57+0.07 

-22.15 

10.24 

10800 

80054 

080916A 

0.6887 

22:25:06.22 

-57:01:22.6 

22.81+0.10 

-19.73 

9.12 

1600 

90062 

080928 

1.6919 

06:20:16.83 

-55:11:58.7 

>24.20 

> -20.24 

< 9.29 

1600 

80153 

081008 

1.967 

18:39:49.87 

-57:25:53.0 

>24.75 

> -19.98 

< 9.18 

5400 

90062 

081029 

3.8479 

23:07:05.36 

-68:09:19.7 

>25.61 

> -20.35 

< 9.39 

10800 

80054 

081109A 

0.9787 

22:03:09.57 

-54:42:40.4 

20.79+0.02 

-22.51 

10.35 

500 

80153 

081118 

2.58 

05:30:22.22 

-43:18:05.6 

25.27+0.41 

-19.98 

9.18 

3600 

80054 

081121 

2.512 

05:57:06.14 

-60:36:09.8 

25.09+0.34 

-20.11 

9.24 

3600 

80054 

081128 

<3.4 

01:23:13.09 

+38:07:38.8 

25.20+0.31 



5400 

90062 

081210 

2.0631 

04:41:56.17 

-11:15:26.8 

22.63+0.04 

-22.20 

10.23 

5400 

90062 

081221 

2.26 

01:03:10.17 

-24:32:52.0 

21.78+0.03 

-23.22 

10.78 

1600 

70036 

081222 

2.77 

01:30:57.57 

-34:05:42.5 

24.48+0.21 

-20.90 

9.61 

3600 

80054 

090313 

3.375 

13:13:36.21 

+08:05:49.8 

>23.54 

> -22.19 

< 10.26 

1600 

80153 

090404 

3.ot? ; | 

15:56:57.49 

+35:30:57.5 

22.14+0.04 

-23.38 

10.85 

1600 

70036 

090417B 

0.345 

13:58:46.61 

+47:01:04.5 

20.94+0.02 

-20.05 

9.44 

500 

70036 

090418A 

1.608 

17:57:15.17 

+33:24:20.9 

23.39+0.11 

-20.95 

9.61 

5400 

90062 

090424 

0.544 

12:38:05.11 

+ 16:50:14.4 

21.35+0.03 

-20.66 

9.62 

1600 

90062 

090516A 

4.109 

09:13:02.61 

-11:51:14.6 

22.92+0.04 

-23.15 

10.71 

10800 

80054 

090519 

3.85 

09:29:07.00 

+00:10:48.9 

>25.50 

> -20.46 

< 9.44 

10800 

80054 

090530 

1.266 

11:57:40.50 

+26:35:37.9 

>24.34 

> -19.51 

< 8.93 

5400 

90062 

090618 

0.54 

19:35:58.48 

+78:21:24.7 

22.63+0.07 

-19.36 

9.04 

1600 

90062 

090709A 

1 q40.5 

19:19:42.63 

+60:43:39.3 

23.58+0.09 

-20.98 

9.63 

4500 

70036 

090715B 

3.00 

16:45:21.61 

+44:50:20.4 

23.45+0.06 

-22.07 

10.18 

10800 

80054 

090812 

2.452 

23:32:48.56 

-10:36:17.2 

>24.80 

> -20.35 

< 9.35 

3600 

80054 

090814A 

0.696? 

15:58:26.39 

+25:37:53.2 

22.98+0.09 

-19.58 

9.05 

1600 

90062 

090926B 

1.24 

03:05:13.94 

-39:00:22.5 

21.42+0.02 

-22.38 

10.28 

5400 

90062 

091018 

0.971 

02:08:44.74 

-57:32:54.1 

22.28+0.04 

-21.01 

9.62 

2700 

90062 

091029 

2.752 

04:00:42.62 

-55:57:20.0 

>24.03 

> -21.34 

< 9.81 

3600 

80054 

091109A 

3.076 

20:37:01.81 

-44:09:30.1 

25.49+0.35 

-20.08 

9.25 

10800 

80054 

091127 

0.490 

02:26:19.87 

-18:57:08.6 

23.21+0.12 

-18.57 

8.66 

1600 

90062 

091208B 

1.0633 

01:57:34.11 

+ 16:53:22.9 

>25.22 

> -18.26 

< 8.28 

2700 

90062 

100305A 


11:13:28.07 

+42:24:14.3 

>25.27 



5400 

90062 

100615A 

1.398 

11:48:49.34 

-19:28:51.5 

23.84+0.14 

-20.21 

9.27 

4500 

70036 

100621A 

0.542 

21:01:13.02 

-51:06:22.4 

21.35+0.03 

-20.65 

9.61 

900 

70036 

100728B 

2.106 

02:56:13.46 

+00:16:52.1 

>24.74 

> -20.13 

< 9.25 

3600 

80054 

100802A 

<3.1 

00:09:52.43 

+47:45:18.5 

25.47+0.43 



5400 

90062 

100814A 

1.44 

01:29:53.55 

-17:59:43.7 

23.35+0.07 

-20.76 

9.52 

5400 

90062 

110205A 

2.22 

10:58:31.20 

+67:31:31.0 

23.80+0.10 

-21.17 

9.72 

11100 

90062 

110709B 

2.09? 

10:58:37.11 

-23:27:16.8 

24.82+0.31 

-20.03 

9.20 

5400 

90062 

120119A 

1.728 

08:00:06.95 

-09:04:53.6 

22.89+0.07 

-21.59 

9.91 

5400 

90062 

120308A 

<3.7 

14:36:20.12 

+79:41:11.9 

24.98+0.27 



5400 

90062 


a Location of IRAC aperture center (J2000). For targets 
with a detected host at any waveband this is the host galaxy 
centroid; for non-detected hosts the best afterglow position 
is given. 

b Measured IRAC aperture magnitude (3.6/rm, AB) and la 
uncertainty. 

c Absolute AB host magnitude at a wavelength of 
A =3.6/mi/(l+2). For targets with photometric redshifts 
we assume the best-fit redshift given in the table. 
d Logarithm of the stellar mass, as derived from the Spitzer- 
measured luminosity. 
e IRAC exposure time in seconds. 

f Spitzer Program ID number. Programs are: 272- PI 
Le Floc’h, “Probing the Dark Side of the Cosmic Evo¬ 
lution Using Dark Gamma-Ray Bursts”. 3653-PI Gar- 
navich, “Gamma-Ray Burst Physics in the Spitzer/Swift 
Era”. 20000 and 30000-PI Berger, “Gotcha! Using Swift 
GRBs to Pinpoint the Highest Redshift Galaxies”. 40598- 
PI Jones, “GRBs as Beacons of Star Formation at High 
Redshifts”. 40599-PI Chary, “Unveiling the Galaxy Coun¬ 
terparts of Damped Lyman-alpha Absorbers using GRB- 
DLAs”. 50562-PI Levan, “The nature of dark gamma- 
ray burst host galaxies”. 70036-PI Perley, “The Host 
Galaxies of Dust-Obscured Gamma-Ray Bursts”. 80054- 
PI Berger, “Probing the z>2 Mass-Metallicity Relationship 
with Gamma-Ray Bursts”. 80153-PI Perley, “Understand¬ 
ing the Environmental Dependence of High-z Dust with GRB 



18 


Per ley et al. 


REFERENCES 


Asplund, M., Grevesse, N., Sauval, A. J., Sz Scott, P. 2009, 
ARA&A, 47, 481 

Blain, A. W., Si Natarajan, P. 2000, MNRAS, 312, L35 
Bloom, J. S., Kulkarni, S. R., Sz Djorgovski, S. G. 2002, AJ, 123, 
1111 

Boissier, S., Salvaterra, R., Le Floc’h, E., et al. 2013, A&A, 557, 
A34 

Bouwens, R. J., Illingworth, G. D., Franx, M., et al. 2009, ApJ, 
705, 936 

Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2014, ApJ, 
793, 115 

Bromberg, O., Nakar, E., Sz Piran, T. 2011, ApJ, 739, L55 
Bruzual, G., Si Chariot, S. 2003, MNRAS, 344, 1000 
Butler, N. R., Bloom, J. S., Sz Poznanski, D. 2010, ApJ, 711, 495 
Cameron, E. 2011, PAS A, 28, 128 

Castellano, M., Fontana, A., Grazian, A., et al. 2012, A&A, 540, 
A39 

Castro Ceron, J. M., Michalowski, M. J., Hjorth, J., et al. 2010, 
ApJ, 721, 1919 

Chabrier, G. 2003, PASP, 115, 763 

Chen, H.-W., Perley, D. A., Pollack, L. K., et al. 2009, ApJ, 691, 
152 

Cucchiara, A., Levan, A. J., Fox, D. B., et al. 2011, ApJ, 736, 7 
Djorgovski, S. G., Frail, D. A., Kulkarni, S. R., et al. 2001, ApJ, 
562, 654 

Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006, ApJ, 644, 813 
Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 
Finkelstein, S. L., Papovich, C., Salmon, B., et al. 2012, ApJ, 756, 
164 

Fryer, C. L., Sz Heger, A. 2005, ApJ, 623, 302 
Fynbo, J. P. U., Prochaska, J. X., Sommer-Larsen, J., Dessauges- 
Zavadsky, M., Sz Mpller, P. 2008, ApJ, 683, 321 
Fynbo, J. P. U., Starling, R. L. C., Ledoux, C., et al. 2006, A&A, 
451, L47 

Fynbo, J. P. U., Jakobsson, P., Prochaska, J. X., et al. 2009, ApJS, 
185, 526 

Galametz, A., Grazian, A., Fontana, A., et al. 2013, ApJS, 206, 10 
Graham, J. F., Sz Fruchter, A. S. 2013, ApJ, 774, 119 
Greiner, J., Fox, D. B., Schady, P., et al. 2015, ApJ, 809, 76 
Hashimoto, T., Perley, D. A., Ohta, K., et al. 2015, ApJ, 806, 250 
Heger, A., Langer, N., Sz Woosley, S. E. 2000, ApJ, 528, 368 
Hirschi, R., Meynet, G., Sz Maeder, A. 2005, ASzA, 443, 581 
Hunt, L. K., Palazzi, E., Michalowski, M. J., et al. 2014, A&A, 565, 
A112 

Izzard, R. G., Ramirez-Ruiz, E., Sz Tout, C. A. 2004, MNRAS, 348, 
1215 

Jakobsson, P., Hjorth, J., Malesani, D., et al. 2012, ApJ, 752, 62 
Kajisawa, M., Ichikawa, T., Tanaka, I., et al. 2011, PASJ, 63, 379 
Kelly, P. L., Filippenko, A. V., Modjaz, M., Sz Kocevski, D. 2014, 
ApJ, 789, 23 

Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35 
Kewley, L. J., Sz Ellison, S. L. 2008, ApJ, 681, 1183 
Kistler, M. D., Yiiksel, H., Beacom, J. F., Sz Stanek, K. Z. 2008, 
ApJ, 673, L119 

Kobulnicky, H. A., Sz Kewley, L. J. 2004, ApJ, 617, 240 
Kocevski, D., West, A. A., Sz Modjaz, M. 2009, ApJ, 702, 377 
Kruhler, T., Greiner, J., Schady, P., et al. 2011, A&A, 534, A108 
Kruhler, T., Malesani, D., Milvang-Jensen, B., et al. 2012, ApJ, 
758, 46 

Kruhler, T., Malesani, D., Fynbo, J. P. U., et al. 2015, A&A, 581, 
A125 

Langer, N., Sz Norman, C. A. 2006, ApJ, 638, L63 
Laskar, T., Berger, E., Sz Chary, R.-R. 2011, ApJ, 739, 1 
Le Floc’h, E., Charmandaris, V., Forrest, W. J., et al. 2006, ApJ, 
642, 636 

Le Floc’h, E., Due, P.-A., Mirabel, I. F., et al. 2003, ASzA, 400, 
499 

Levesque, E. M. 2014, PASP, 126, 1 

Levesque, E. M., Kewley, L. J., Graham, J. F., Sz Fruchter, A. S. 
2010, ApJ, 712, L26 

Lien, A., Sakamoto, T., Gehrels, N., et al. 2014, ApJ, 783, 24 
MacFadyen, A. I., Si Woosley, S. E. 1999, ApJ, 524, 262 
Madau, P., Sz Dickinson, M. 2014, ARA&A, 52, 415 


Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463 
Modjaz, M., Kewley, L., Kirshner, R. P., et al. 2008, AJ, 135, 1136 
Muzzin, A., Marchesini, D., Stefanon, M., et al. 2013, ApJS, 206, 
8 

Niino, Y. 2011, MNRAS, 417, 567 

Oesch, P. A., Labbe, I., Bouwens, R. J., et al. 2013, ApJ, 772, 136 
Peng, C. Y., Ho, L. C., Impey, C. D., Si Rix, H.-W. 2002, AJ, 124, 
266 

Perley, D. A., Levan, A. J., Tanvir, N. R., et al. 2013, ApJ, 778, 
128 

Perley, D. A., Perley, R. A., Hjorth, J., et al. 2015a, ApJ, 801, 102 
Perley, D. A., Kruhler, T., Schulze, S., et al. 2016, ApJ, 817, 7 
Pettini, M., Sz Pagel, B. E. J. 2004, MNRAS, 348, L59 
Podsiadlowski, P., Ivanova, N., Justham, S., Sz Rappaport, S. 2010, 
MNRAS, 406, 840 

Porciani, C., Si Madau, P. 2001, ApJ, 548, 522 
Prochaska, J. X., Chen, H.-W., Dessauges-Zavadsky, M., Sz Bloom, 
J. S. 2007, ApJ, 666, 267 

Prochaska, J. X., Sheffer, Y., Perley, D. A., et al. 2009, ApJ, 691, 
L27 

Ramirez-Ruiz, E., Trentham, N., Sz Blain, A. W. 2002, MNRAS, 
329, 465 

Robertson, B. E., Si Ellis, R. S. 2012, ApJ, 744, 95 
Robertson, B. E., Ellis, R. S., Furlanetto, S. R., Sz Dunlop, J. S. 
2015, ApJ, 802, L19 

Rossi, A., Klose, S., Ferrero, P., et al. 2012, ASzA, 545, A77 
Salvaterra, R., Della Valle, M., Campana, S., et al. 2009, Nature, 
461, 1258 

Sanders, R. L., Shapley, A. E., Kriek, M., et al. 2015, ApJ, 799, 
138 

Santini, P., Ferguson, H. C., Fontana, A., et al. 2015, ApJ, 801, 97 
Savaglio, S. 2012, Astronomische Nachrichten, 333, 480 
Savaglio, S., Glazebrook, K., Si Le Borgne, D. 2009, ApJ, 691, 182 
Savaglio, S., Glazebrook, K., Le Borgne, D., et al. 2005, ApJ, 635, 
260 

Schady, P., Savaglio, S., Muller, T., et al. 2014, A&A, 570, A52 
Schulze, S., Chapman, R., Hjorth, J., et al. 2015, ApJ, 808, 73 
Steidel, C. C., Rudie, G. C., Strom, A. L., et al. 2014, ApJ, 795, 
165 

Svensson, K. M., Levan, A. J., Tanvir, N. R., Fruchter, A. S., Sz 
Strolger, L.-G. 2010, MNRAS, 405, 57 
Svensson, K. M., Levan, A. J., Tanvir, N. R., et al. 2012, MNRAS, 
421, 25 

Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, Nature, 461, 
1254 

Tanvir, N. R., Levan, A. J., Fruchter, A. S., et al. 2012, ApJ, 754, 
46 

Totani, T. 1997, ApJ, 486, L71 

Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 
613, 898 

Trenti, M., Perna, R., Sz Jimenez, R. 2015, ApJ, 802, 103 
Trenti, M., Perna, R., Levesque, E. M., Shull, J. M., Sz Stocke, J. T. 
2012, ApJ, 749, L38 

van den Heuvel, E. P. J., Sz Portegies Zwart, S. F. 2013, ApJ, 779, 
114 

Vergani, S. D., Salvaterra, R., Japelj, J., et al. 2015, A&A, 581, 
A102 

Wanderman, D., Si Piran, T. 2010, MNRAS, 406, 1944 
Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 
1 

Wijers, R. A. M. J., Bloom, J. S., Bagla, J. S., Sz Natarajan, P. 
1998, MNRAS, 294, L13 

Wolf, C., Si Podsiadlowski, P. 2007, MNRAS, 375, 1049 
Woosley, S. E., Si Bloom, J. S. 2006, ARA&A, 44, 507 
Woosley, S. E., Sz Heger, A. 2006, ApJ, 637, 914 
Yoon, S.-C., Si Langer, N. 2005, ASiA, 443, 643 
Yoon, S.-C., Langer, N., Sz Norman, C. 2006, ASiA, 460, 199 
Zahid, H. J., Dima, G. I., Kudritzki, R.-P., et al. 2014, ApJ, 791, 
130 

Zemcov, M., Smidt, J., Arai, T., et al. 2014, Science, 346, 732 



